Fulde-Ferrell-Larkin-Ovchinnikov State and Field-Induced 
Superconductivity in an Organic Superconductor 



m 
o 
o 

(N 



(N 
O 

o 

Of 



^3 

o 
o 



> 

On 
O 

m 
O 
(N 
O 

-i— > 
c3 



i 

c 

o 
o 



X 



Hiroshi Shimahara 

Department of Quantum Matter Science, ADSM, Hiroshima University, Higashi- Hiroshima 739-8530, Japan 

(Received 7 March 2002) 

An experimental phase diagram of the field-induced superconductivity in the A-(BETS)2FeCl4 
compound is theoretically reproduced by a combination of the Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) state and the Jaccarino-Peter mechanism. Semi-quantitative agreement is obtained in the 
presence of a very weak subdominant triplet pairing interaction between electrons with antiparallel 
spins, in addition to a dominant singlet pairing interaction. An order-parameter mixing effect 
enhances the FFLO state and plays an essential role in the semiquantitative agreement. As a 
possible origin of the small triplet component, pairing interactions mediated by spin fluctuations 
are proposed, ft is argued that, in this mechanism of the pairing interactions, equal spin pairing 
interactions are much weaker than the antiparallel spin pairing interactions in the strong magnetic 
field, which is consistent with the lack of observation of equal spin state. 



PACS numbers: 74.80.Dm,74.80.-g 

Recently, a field-induced superconductivity (FISC) 
was observed at high magnetic fields in an or- 
ganic superconductor A-(BETS) 2 FeCl 4 , (BETS = bis- 
(ethylenedithio)tetraselenafulvalene) [1]. The lower crit- 
ical field H l ° w (T) of the FISC is upwards convex [2,3] 
as a function of temperature T. When we consider 
the Jaccarino-Peter mechanism [4] as the origin of the 
FISC [2,5,6], the lower critical field corresponds to the 
upper critical field of a typical superconducting phase. 
Therefore, the upwards convex H l ° w (T) is not a com- 
mon characteristic of the upper critical field. Some au- 
thors argued [2,3,5] that this upturn behavior can be ex- 
plained by the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
state [7-12]. 

The FFLO state has been discussed in some other 
organic superconductors, because some of them satisfy 
the necessary conditions for the FFLO state to occur: 
(1) Those samples can be clean type II superconductors 
due to their narrow electron bands; (2) The orbital pair- 
breaking effect is strongly suppressed for parallel mag- 
netic fields due to the low dimensionality. The present 
compound A-(BETS)2FeCi4 also exhibits these features. 

The upturn due to the FFLO state in the cylindrically 
symmetric system was obtained by Bulaevskii [9] and Aoi 
et al. [10]. The FFLO state in two-dimensional supercon- 
ductors was reexamined and discussed in connection with 
exotic superconductors, such as organics and cuprates, 
by some authors [11-19]. The upturn is a common be- 
havior of the FFLO state in quasi-two-dimensional sys- 
tems [13-19]. This behavior is due to a Fermi surface 
effect [11,13,14] analogous to Fermi surface nesting for 
the spin-density and charge-density waves, which we call 
a nesting effect for the FFLO state. 

On the other hand, an upturn of the upper critical 
field (d 2 H c /dT 2 > 0) and a first-order transition below 
the upper critical field have been observed in k-(BEDT- 
TTF) 2 Cu(NCS) 2 [20-23]. The FFLO state has also been 
discussed in this compound. Manalo and Klein [24] have 
successfully confirmed the agreement between the exper- 



imental [20-22] and theoretical results [17]. 

A recent experiment on thermal conductivity by 
Tanatar et al. [25] supports the existence of the FFLO 
state in A-(BETS)2GaCLi. They observed that the up- 
per critical field is consistent with the Pauli paramag- 
netic limit in dirty samples, while it shows no satura- 
tion down to 0.3 K in clean samples. This experimental 
result supports the FFLO state also in the FISC in A- 
(BETS)2FeCl4, since these two compounds are similar, 
except A-(BETS)2FeCi4 has the magnetic anions. 

In A-(BETS) 2 FeCl 4 , localized spins of S = 5/2 on the 
(FeCU) -1 anions create an exchange field on the con- 
ducting layer of the two-dimensional network of BETS 
molecules. As many authors pointed out [1,2,5,6], the 
exchange field shifts the Zeeman energy effectively, and 
thus also the superconductivity area in the B-T phase 
diagram. This compensation mechanism is called the 
Jaccarino-Peter mechanism, which was originally pro- 
posed for the FISC in ferromagnetic metals. The present 
compound is an antiferromagnetic insulator at zero field, 
but this difference does not affect the mechanism at high 
magnetic fields, where the localized magnetic moments 
are saturated [2,5,6]. 

Balicas et al. have compared the experimental phase 
diagram of A-(BETS)2FeCU with the theoretical one 
with the FFLO state, taking into account the shift of the 
FISC area by the Jaccarino-Peter mechanism [2]. The 
main features of the FISC have been reproduced qual- 
itatively by a combination of the FFLO state and the 
Jaccarino-Peter mechanism. 

However, from a quantitative viewpoint, the theoret- 
ical width of the magnetic field region of the FISC is 
much narrower than the experimental one. In the ex- 
perimental phase diagram, the FISC area has a center 
at about 33 T and the lower limit at about 17 T [3], 
which means that the absolute value of the critical field 
is equal to about 16T in the absence of the exchange 
field. When the orbital pair breaking effect is ignored, 



the maximum transition temperature T ( 



(m) 



4.2 K gives 
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the Pauli paramagnetic limit 1.856[T/K] x T< 



(m) 



and the s-wave FFLO critical field 11 T for the cylin- 
drical Fermi surface. Therefore, the experimental value 
16 T is much larger than the estimates using the simple 
theory. 

Therefore, we need an explanation for this difference. 
There are some mechanisms which enhance the FFLO 
critical field, such as the effects of the anisotropies of the 
Fermi surface and the gap function [13-18]. In particular, 
by a nesting effect, the critical field can be enhanced up to 
a value several times as large as the Pauli paramagnetic 
limit [13,14]. However, the increase due to this mecha- 
nism occurs at low temperatures [15], and the tricritical 
point does not change. In the case of A-(BETS) 2 FeCl4, 
it seems that the shape of the FISC area in the experi- 
mental phase diagram is not reproduced very well solely 
by these mechanisms. 

Another mechanism is an order-parameter mixing ef- 
fect. In the presence of a subdominant triplet pairing in- 
teraction between electrons with antiparallcl spins, even 
if it is very weak, the upper critical field is largely en- 
hanced. This effect was theoretically predicted in s-wave 
superconductors with a spherical symmetric Fermi sur- 
face [26] , and also in s- and <i-wave superconductors with 
quasi-two-dimensional Fermi surfaces [19]. In the latter 
study, this effect was discussed in connection with the 
organic superconductors. Subdominant triplet pairing 
components exist in the pairing interactions mediated by 
the spin fluctuations [27,28], and also in those mediated 
by phonons [29]. In the latter, the triplet pairing com- 
ponent is enhanced due to the weakness of the screening 
effect in the low dimensions [30]. 

In this paper, we reproduce a phase diagram that 
agrees with the experimental data semiquantitatively, by 
taking into account the order-parameter mixing effect in 
the FFLO state with the Jaccarino-Peter mechanism. 

Before we analyze of the FISC, let us briefly discuss 
the exchange field created by the (FeCU) 1 anions. In 
our previous paper [5], we proposed a mechanism which 
enhances the upper critical field in antiferromagnetic su- 
perconductors. In order to illustrate the mechanism, we 
used a generalized Kondo lattice model of the localized 
spins and mobile electrons with a Kondo coupling Jk be- 
tween them and an antiferromagnetic exchange coupling 
J between the localized spins. As a special case, the 
model describes an aspect of the A-(BETS)2FeCLt com- 
pound, for a half-filled electron band and large ratios of 
Jk/zJ, where z is the number of the nearest neighbor 
sites. In the application to the A-(BETS)2FcCl4 com- 
pound, J should be replaced with some indirect inter- 
action via conducting layers [6,31], which must be much 
smaller than Jk. Therefore, it is expected that the ratio 
Jk/zJ is much smaller than 1. The RKKY exchange in- 
teraction and the reconstructed Fermi surface have been 
calculated by Brossard et al. [31]. 

In the calculation based on the generalized Kondo 
lattice model [5], the boundary between the antifcrro- 



7.8 T magnetic insulating phase and the paramagnetic metal- 
lic phase is Bmi = HqHmi = ^qzJS/\^ c \, where /xo is the 
magnetic permeability and /j, c is the electron magnetic 
moment. On the other hand, the compensation effect 
shifts the center of the FISC area by 

Scent = A*0 H con t = fl -r^-r. (1) 
lMe| 

Therefore, from the experimental data [1], -Bmi ~ 10. 5T 
and i? ccn t ~ 33T, we obtain the ratio 



Jk 

zJ 



Bmi 



3.1 > 1 , 



(2) 



which is consistent with the above consideration. 

From the values -B C ent = 33 T, Bmi = 10.5 T and 
S = 5/2, we obtain J K /k B « 8.87 K and z J/k B = 2.82 K. 
The mean field value of the antiferromagnetic transition 
temperature Tf F F is equal to zJS(S+l)/4 « 8.1 K, which 
agrees very well with the experimental value 8.5 K. 

Now, let us examine the FISC. We examine the case 
in which a dominant singlet pairing interaction and a 
very weak subdominant triplet pairing interaction coexist 
between electrons with antiparallel spins. The pairing 
interaction and the gap function are expanded as 



V(k,k') = -5> Q7 a(k)7a(k'), 
A(k) - ]TA Q7Q (k), 



(3) 
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with symmetry factors 7«(k). In cylindrically symmet- 
ric systems, they are defined by, for example, 7 s (k) oc 1, 
7 Px (k) oc k x , ld x 2_ 2 (k) oc k 2 x — k y , and so on. We nor- 
malize 7 Q (k) by 

f -2lT 

^[7-(v)] a = l, (4) 

where J a (f) — [7a(k)]iki=fe P and ip is the angle between k 
and k x -xas. We assume a cylindrically symmetric Fermi 
surface for simplicity, and s-wave and p-wave pairing in- 
teractions as the dominant and subdominant pairing in- 
teractions, respectively [34]. Later, we will discuss the 
effects of the anisotropies of the Fermi surface and the 
gap function. 

In the weak coupling limit, the gap equations are writ- 
ten as 

T 



A a log 
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where we define 
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" y (cosh 2 y + sinh 2 4r) 



= 2eWe- 1 /^(°), 



(6) 
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with ( = h [qcos((p — tpq) + l], q = v F q/2h, h = \fj, e \(\H\- 
H C ent): N(0) the density of states at the Fermi energy, 
and ifq the angle between the center-of-mass momentum 
q of the FFLO state and the k x -&xis. 

For the symmetry, we can assume the direction of 
q along the fc^-axis without loosing generality. Then, 
the Py-wave component A p k v is not mixed with the s- 
wave component, while the p^-wave component A Pn .k x is 
mixed with it. Therefore, the transition temperature for 
a given q is calculated by T c (/i,q) = T c ( ° ) e _A (' l ' q ) with 
the smallest eigenvalue X(h, q) of the 2x2 matrix 



M, 



M x . 



M ps M pp + G p 



where we define 



,(o) 



G p = log ■ 



r (o) 



g p N(0) g s N(0)- 



(7) 



(8) 



Here, and Tc P are the pure s- and p-wave transition 
temperatures. In the presence of the exchange field, 
is equal to the maximum transition temperature ri m ' of 
the FISC. The final result of the transition temperature 
is obtained by maximizing T c (h, q) with respect to q. 

The results are depicted in Fig. 1. The FFLO crit- 
ical fields are largely enhanced even for a very small 
value of the coupling constant of the coexisting triplet 
pairing interaction. Figure 1(a) shows the result for 
T c (m) = 4.2 K and = 0.084 K, which corresponds to 

Tc ( pV T c (m) = 1/50 and G p w 3.912. It is found that the 
agreement between theoretical and experimental results 
is significantly improved by taking into account the weak 
triplet pairing interaction (solid curve), in comparison to 
the result in the absence of the triplet pairing interactions 
(dotted curve) . Taking into account the ambiguity in de- 
termination of the experimental values of the transition 
temperatures [1-3], it can be concluded that the over- 
all behavior of the transition curve is semiquantitatively 
reproduced by the FFLO state. 

Figure 1(b) is the result for T c (m) = 4.5 K and T$ = 



,(o) 



0.045 K, which corresponds to T^> /T c 



(0) ,™(m) 



1/100 and 



G p « 4.605. Except near the top of the peak, the theo- 
retical curve agrees very well with the experimental data. 

In the above calculations, the presence of the weak 
triplet pairing interactions between electrons with an- 
tiparallel spins is assumed. When the system is isotropic 
in the spin space, triplet pairing interactions between 
electrons with parallel spins are of the same magnitude 
at zero field. In this case, one might expect transitions 
to an equal spin pairing state at magnetic fields where 
the FISC of antiparallel spin pairing is suppressed. For 
example, in the cases shown in Figs. 1(a) and 1(b), the 
transition temperatures of the equal spin state are equal 
to 84 mK and 45 mK, respectively. However, such tran- 
sitions have not been observed in the experiments. 
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FIG. 1. Comparison between the theoretical and experi- 
mental phase diagrams in A-(BETS)2FeCl4. We set the max- 
imum value of the superconducting transition temperature 



(a) T t 



(m) 



4.2 K and T$ = 0.084 K, and (b) T< 



(m) 



4.5 K 



and Tcp^ = 0.045 K. The strength of the exchange field is 
assumed so that the center of the peak B C cnt is equal to 33 T. 
The solid, dotted, and dashed curves show boundaries of the 
FISC, in the presence of the triplet pairing interaction, in 
the absence of it, and when the FFLO state is excluded, re- 
spectively. The solid circles and the solid triangles show the 
experimental data of the superconducting and antiferromag- 
netic transition points, respectively, by Uji et al. in ref. [3]. 
The dot-dashed curve is a guide for the eyes. 

This discrepancy can be explained if the origin of 
the triplet pairing interactions is the exchange of the 
spin fluctuations enhanced by on-site Coulomb repul- 
sion [27,32]. For example, in a random phase approxi- 
mation (RPA), the equal spin pairing interactions orig- 
inate from the fluctuations as shown in the diagrams in 
Fig. 2(a) [27]. Therefore, they correspond to the fluc- 



3 



tuations of S z components. On the other hand, the 
antiparallel spin pairing interactions originate from the 
fluctuations as shown in the diagrams in Figs. 2(b) and 
2(c) [27,32] in the RPA. They correspond to the fluctua- 
tions of S z components and S± = S x ± iS y , respectively. 
In the magnetic field in the ^-direction, since the spin 
moment (S z ) is finite, the spin fluctuations of z compo- 
nents are much smaller than those of x and y components. 
Therefore, the equal spin pairing interactions are much 
weaker than the antiparallel spin pairing interactions in 
the magnetic field. 



(a) iQ + Qt + - 



upwards convex (d 2 H c (T)/dT 2 < 0) can be explained by 
the combination of the FFLO state and the Jaccarino- 
Peter mechanism, together with the order-parameter 
mixing effect. For the semiquantitative reproduction of 
the phase diagram in the present theory, the existence of 
the very weak subdominant triplet pairing interaction is 
essential. The origin of it can be attributed to spin fluc- 
tuations in the magnetic field. The phase diagrams of 
A-(BETS)2Gai_ x Fe x CLi are also expected to be under- 
stood in the present mechanism. The experimental evi- 
dence of the existence or inexistence of the FFLO state 
can be obtained by direct observations of the structure 
of the gap function [33,35]. 

The author would like to thank Dr. Uji for useful dis- 
cussions and the experimental data, and Dr. Cepas for 
useful discussions. 
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FIG. 2. Series of diagrams of fluctuations which contribute 
to the pairing interactions. The short dashed line denotes the 
on-site Coulomb interaction U. 



Another possible explanation is the orbital pair- 
breaking effect. The transition temperature of the equal 
spin pairing state is only 45 mK in the case shown in 
Fig. 1(b). If we take into acount the orbital pair-breaking 
effect, the equal spin pairing state with such a small 
transition temperature may be suppressed completely at 
higher magnetic fields where the system is metallic. 

In this paper, we have assumed the s-wave pairing in- 
teraction as the dominant singlet interaction, and the 
cylindrical Fermi surface for simplicity [34]. For the d- 
wave pairing, the curve of the critical field exhibits a 
kink at a low temperature [16-18]. However, the kink 
disappears when the Fermi surface is anisotropic to some 
extent [15]. Furthermore, the anisotropy of the Fermi 
surface enhances the critical field [13,14] for a nesting ef- 
fect. However, since the enhancement occurs especially 
at low temperatures T < 0.2T C [15], the mixing effect re- 
mains necessary for the semiquantitative reproduction of 
the experimental phase diagram in an entire temperature 
region including T > 0.2T C . 

In conclusion, the phase diagram of A-(BETS)2FeCl4 
including the curve of the lower limit of the FISC which is 
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